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The 77-3./V-interaction is studied within the four-body scattering theory adopting purely separable 
forms for the two- and three-body subamplitudes, limiting the basic two-body interactions to s- waves 
only. The corresponding separable approximation for the integral kernels is obtained by using the 
Hilbert-Schmidt procedure. Results are presented for the ??- 3 H scattering amplitude and for the 
total elastic cross section for energies below the triton break-up threshold. 
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I. INTRODUCTION 

In the last ten years the interaction of 77-mesons with few-nucleon systems has attracted considerable interest, and 
quite a large amount of research has been carried out both on the experimental as well as on the theoretical side. The 
major aim of this activity is to obtain a well-defined conceptual picture of the low-energy ry-nuclear interaction which 
then may serve to clarify the foundation of the more fundamental r)N problem. 

Of special interest is the rj-3N -system, for which exactly soluble models are available. At the same time it encloses 
a much wider range of phenomena than the more simple r/NN case. In this context we would also like to mention 
several precise experimental investigations of r\- production on three-body nuclei performed during the last decade . 
Previous theoretical studies in this area were restricted to various approximations in treating the few-particle aspects 
Cn| ' of this problem. As an earlier nonperturbative calculation, we would like to refer to the optical model approach of [|J, 
whose purpose was to estimate qualitatively the possible formation of bound ?7-meson states with the lightest nuclei. 
More refined calculations were presented in where the authors were able to sum the multiple scattering series 
for the ry- 3 H amplitude, including several important corrections to the trivial optical limit. Another result reported 
in H were obtained within the so-called finite-rank-approximation (FRA) which has recently been applied also to the 
\ ^-production from three-body nuclei . The crucial point of this model is the neglect of target excitations during 
the interaction with the ^-meson. Clearly this assumption allows one to avoid the complications associated with the 
direct solution of the four-body dynamical equations. Concerning the present study, we were guided by the idea 
that the approximations used in many-body physics, such as the optical model or adiabatic treatment of the target, 
may fail when few-particle systems are studied. Especially, this seems to be true for small kinetic energies, where 
• i-H , the unitarity conditions of the scattering matrix, nucleon recoil and other effects become significant. Their neglect, 
tolerated in various many-body approaches, may affect drastically the quality of the few-body results. Therefore, the 
present work is intended to solve the 77-37V-problem without making any such not well controlled approximations. 

The basis of our calculation is the four-body formalism in momentum space. Although a variety of methods for 
solving the n-particle problem has been proposed in the literature, the Faddeev-Yakubovsky theory and the one 

of Alt, Grassberger and Sandhas (AGS) |ll] : l^] are most convenient and preferred for practical applications. Adopting 
the separable representation for the driving two-body potentials as well as for the subamplitudes appearing in the 
(1+3)- and (2+ 2)-p artitions of the four-body system, both approaches lead to the same set of effective two-body 
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equations [[11 13 14[] . As is well known, the separable approximation of the integral kernels permits one to represent 
the dynamical equations in terms of particle exchange diagrams. Due to this tractability and its relatively simple 
numerical realization, this method has received wide acceptance in few-body physics. Thus at present, a feasible 
formalism of four-particle theory has been extensively developed, despite the much more complex structure of the 
corresponding equations compared to the three-body case. After the work of Tjon fl5| impressive results have been 
obtained in recent years for the four- nucleon low-energy interaction (see e.g. |l(| and references therein), as well as for 
pion absorption on three- nucleon systems Jl9|-pl|. With respect to other techniques we would like to refer to recent 
work in lO and El. 
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The paper is organized as follows. In the next section we outline briefly the formal aspects concerning the application 
of the four-body formalism to the tj-3N system within the quasiparticle approach. Besides the basic equations we 
introduce here the two- and three-body ingredients of the model. Our results for ?y- 3 H elastic scattering are presented 



in Sect. Ill, and conclusions are drawn in the last section. Details are given in two appendices. 



II. APPLICATION OF THE FOUR-BODY FORMALISM TO THE T]-3N SYSTEM 

The separable method is well known to allow one to reduce the n-body problem to an in general simpler (n — 1)- 
body case, where one of the constituents appears as a quasiparticle, i.e., a two-body bound, virtual or resonance 
state. In particular, the three-body equations for the (3— »3) transition amplitudes are exactly reduced to effective 
two-particle equations of Lippmann-Schwinger type. Their kernels contain the off-energy-shell (2— >2) amplitudes 
for all two-body subsystems. In an analogous manner, the four-body scattering kernels can be expressed in terms 
of subamplitudes stemming from the decomposition of the four-body system into the partitions (1+3) and (2+2). 
Therefore, approximating these subamplitudes once more by a separable ansatz we can again reduce the four-body 
problem to an effective quasiparticle two-body one. The formal details of this two-step reduction scheme may be 
found in Refs. 0|J|,[[|. Here we restrict ourselves to a brief description of the resulting quasi-two-body equations 
as applied to the t]-3N system. 

A. The four-body r/SN equations 

Since our formalism does not include Coulomb forces, we will consider for definiteness the ?7- 3 H interaction through- 
out this paper. All appropriate states are assumed to be properly antisymmetrized with respect to the nucleons. The 
corresponding antisymmetrization procedure is outlined in Appendix A. Then we are led to the following three chan- 
nels, corresponding to three possible two-quasiparticle partitions of the rj-3N system 

(1): 7/ + (3JV) , (2): N + (rjNN) , (3): (Nrj) + (NN) . (1) 

We only need the amplitudes connecting the initial asymptotic state, consisting of the 3N bound state ( 3 H) and a 
free 77- meson, with all three channels listed in ([!]). They obey a set of three coupled integral equations (see Appendix 
A), whose structure is represented by the following matrix equation 

Z12 Z13 

Z21 z 22 z 23 I j e 2 ) I x 2 I . (2) 

Z31 Z32 

Here the index a=l,2,3 stands for the channel (a) from (|l|). The amplitude X a describes the transition (1) — > (a). 
The effective potentials^] Z a p are expressed in terms of the form factors, generated by the separable representation of 
the subamplitudes appearing in the channels (]!]). 

Because here we consider only low energy scattering, we take into account only the dominant s-wave part of the 
interaction in the two-body subsystems and thus only s-waves in the three- and four-particle states. Therefore, the 
matrix elements are diagonal with respect to the total spin S. For the elastic 7?- 3 H scattering we need to consider only 
those states where the spins and isospins of all nucleons are coupled to S = T = 1/2. In all expressions to follow we 
drop the index L — 0. The explicit analytical form of the potentials Z a p, taking into account the spin-isospin degrees 
of freedom, is given in the Appendix B. In a more detailed notation the system (^]) reads 

XLZUP,P\E) = (1 - S al )Z^ n ,(p,p',E) 

o OO 

-(»"*') (II I T?\ P " 2 d P 
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[3—1 n" s" 



(3) 



1 Following the work |HJ we explore the formal analogy with the Lippmann-Schwinger equation and use for the driving terms 
Z a p the suggestive term "potential" . 
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The index s — (0, 1) in the above equations corresponds to the total spin of the given NN pair. Clearly, due to the 
pseudoscalar-isoscalar nature of the 77-meson, the value of s fixes uniquely the spin structure of the overall 4-body 
state with total spin S—l/2. In view of the limitation of the two-body interaction to the dominant s-wave part, the 
isospin t of a NN pair is fixed by its spin s through the condition s + 1 = 1. The subenergies £p in (|^) are defined as 

£,=E-^-,(P = W), (4) 

with reduced masses 

ZM N m v M N (2M N + m n ) 2M N (M N + m v ) 
Mi = , 1VI2 — , IVlz — • \P) 

3-Mjv + irir) 3Mjv + m ri 3Mn + rn. n 

A graphical representation of the system (||) is shown in Figs. |l| and |[ 

The structure of (Q) allows one to eliminate the channel (1) yielding an equivalent set of only two coupled equations 
for the amplitudes X2 and X3. In detail, one has for a € {2, 3} 



00 

R — O 1 a" „ 



p" 2 dp 



13=2,3 n"s" 

where the new effective potentials are given by (a, (3 € {2, 3}) 

Z l :H n ,(p,p',E) = (l - Ml - S Q2 ))zl%l n ,(p,p\E) 



2tt 



2 



(6) 



(7) 



n"s" 

After solving the system (||), the amplitude X\ is obtained from 

oc 

x^,(p, P ',e)= EE/ ^^^m^M^xffi^m^- (8) 



13=2,3 ri's" 



The set of equations (^|) is more suitable for the numerical solution than ([}]), since in the former case the integration 
over the triton pole in the propagator 0i may be carried out independently from the procedure of solving the system 
(^J) itself. Thus the kernels in ([3]) are smooth functions of the integration variable, and the equations may be solved 
by direct matrix inversion. We recall that below the triton break-up threshold we are dealing with only nonsingular 
potentials Z a $. 



B. The subamplitudes 

The key ingredient of the quasiparticle method pl| , p^4[ |, leading to the equations of the type (^), is the separable 
representation of the off-shell scattering amplitudes for the two- and three-body subsystems, appearing in the (2+2)- 
and (l+3)-partitions. In our case, the two types of two-body subsystems involved are NN and r]N, denoted in the 
following by "d" and "jV*", respectively. For the corresponding scattering matrices we adopt the simplest rank-one 
separable form. In detail, we use for the NN interaction 

t l ; ) ( q ,q',z)=g^( q )r^(z) 9 l ;\q'), (9) 

with 

2Mm L 47r^ J zMn — q - 1 


Here the upper index (s) stands for the singlet (s = 0) and triplet (s — 1) NN states. For the form factors (q) we 
use the simplest parameterization 
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Analogously, we choose in the rjN channel 

t N .(q,q',z) = gl,(q)r N ,(z)gl,(q') , (12) 



with 



where 



-1 

t n * (z) = - z + M N + m n - M - %ri(z) - ^tt(z) , (13) 



J z + m„ — m — 77"^ h ze Z7r 

and 

p f \ p P"p M N m p 

The parameters appearing in ( |TT| ) , ( |l3| ) and ([[5]) are summarized in Table H. Those for the TViV interaction were taken 
from the low-energy NN fit presented in p2| , and the ryiV-parameters are chosen as to reproduce the scattering length 
a i)N = (0.75 + i 0.27) fm which agrees with the most recent results p3| , p4| . However, our parameterization gives a 
different value for the effective range parameter ro = (1.95 + «0.07) fm compared to ro = (1.5 + £0.24) fm of [ pi| . 
Unfortunately, we are not able to fit both values a n N and ro simultaneously, which is of course the price one has to 
pay for using the simplest separable ansatz . When comparing our predictions with those of we use also other 
sets of ryiV-parameters chosen in such a way that they lead in each case to the corresponding value of a^jv used there 
(see Table 

Turning now to the general scheme, we have to introduce also the separable representation for the three-body 
subamplitudes, which will then serve as a necessary input for the four-body calculation. For this purpose we apply 
the Hilbert-Schmidt expansion. The main formal aspects of the procedure can be found e.g. in p5] , |26| j. 

The three-nucleon s-wave doublet amplitudes Ui ;ss >, appearing in the channel (1) obey the equation (see e.g. |27| ) 

00 

U Us Aq,q',£) = V hs Aq,q\£)+ £ / V 1;ss „(q, q" , £ ) rj' \£ - U 1;s „ s ,(q", q', £ )^- , s, s> = 0, 1 . (16) 

s"=0,l{ N 

The effective potentials are defined in terms of the form factors gd(k) of the two-nucleon amplitude (j^) as 

v , , -x a ss , y g ^(ig , +Ki)gd'' ) (ig+ , n , n7) 

Vi; SS >(q,q ,£) = — J q2 ql2 „, d(q-q), s,s =0,1, (17) 

_l M N Mff M N 

where £ denotes the total cm. kinetic energy in the three nucleon system, and the matrix of the spin-isospin coefficients 
is given by 

A =U1)- 

The Hilbert-Schmidt expansion for the driving term Vi- ss i reads 

I/ 1;ss ,( g ,g',£) =-J2K(£)utf(q,£)u£\q',£), (19) 

n 

where the functions Un\q,£) are taken as the eigenfunctions of the kernel of equation ( |l6| ) with the eigenvalues A„, 
i.e. 
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<4 s) (^) = X7?T £ f V UsAQ,q\£)r!?\£-§ r )ui*'X q ',£)'^, , = 0,1. (20) 

' s'=0,1q W 



They are normalized according to 

oo 

s =o,i £ 

The separable form of the amplitude can easily be found 



(21) 



U l ., ss ,{q, q ',£)=Y d ^\q,e)Q 1 , n {E)u^\q',£), Q Un (£) = - X " { ^_ - . (22) 

The amplitudes for the ijNN scattering, related to the channel (2) and denoted in the following by U^ij with 
i,j G {d,N*}, are coupled into two independent sets corresponding to two possible spin-isospin r/NN s-wave states 
(5; T) — (0; 1) and (1; 0), denoted by an upper index (s) with s e {0, 1}, each of the form 

oo 



2:. 



(q,q',£) = J V^ N 4q, q '\£)r N {£ - ^)ui%, d (q'', q ',£) q -^ , (23) 



tfgWg,^) = vg. d (<z,<z',£) + I [V 2 % d (q, q ",£)4 s \£ - ^-)u { 2 %(q", q',£ ) 



//2 

+^S. w »(9,9",f)r J v,(f-^-)l7^, d (g",g',£) 



2/ijV* 

The reduced masses appearing in are defined by 



q " 2 d q " 
2vr 2 



„ _ M N (M N + mr,) ~ = 2M N m v 
~ 2M N + m n ' ^ d ~ 2Mjv + m„ ' 



The corresponding effective potentials are 



V 2 \J N . {q,q,£) = -j=J q2 g w d(q ■ q ) , 

_1 Mat 2fi v N Mat 



^, d (<z, g', f ) = V^V (<?', 9, f ) , (25) 

„w ( , i + M»(ig , +^gi)g^(ig+^g , i) J ,. ^ 

^V(9.9.^) = o / - J ^ d{q-q) 



£ - -2— 



>r<i 



where £ is the total cm. kinetic energy of the rjNN system. The properly antisymmetrized amplitudes for rjNN 
scattering are expressed in terms of U 2 .^ as (see I28|j29(| ) 

^2 ;dd = ^2;eL ) ^2;AT*d = ^^^2;Af*d > S = 0, 1 . (26) 



Analogously to < \22\j we have (i, j € {d, N*}) 

^ ; (5)-i 

where the form factors obey the homogeneous equation 
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oo t 



Here, of course, one has rjy» = t^. (s = 0, 1), and the cigenfunctions are normalized as 

oo 

E /^<*^(«-0©(..«>g--*~- <»> 

i=d,N' J Q ^ % 

In the actual calculation we have neglected any ttNN states. Their inclusion would imply an increase in the number of 
channels in the final equations as well as adoption of relativistic kinematics which would lead to much more complex 
formalism. On the other hand due to its small mass, the pion is expected to give only minor corrections to low-energy 
•q- nucleus scattering |p9 30 1. 

Apart from the genuine three particle scattering amplitudes, we also need as input the effective amplitudes, denoted 
here by U^J-, which describe two independent pairs of interacting particles in the channel (3). The corresponding 
equations read in our case 

oo 

/ II2 > n"^rln" 
V 3 %, (q, q",£) r N {£ - ±-) vW.itf', £) i^- , (30) 

o 

/ II2 . Ill J II 

V 3 % d (l> 9", £) r [ d s) (£ - fj Uj&tf, </, £) q -^~ , 



where £ is the sum of the internal energies in the NN and r/N subsystems. In the expressions (|30j) the notations 
Vd = H-qN and vpj* — Mjv/2 are used. The effective potentials are 

r( S ) („ „' P\ - T/( S ) 



v {s) 



V 3 % N ,(q,q J ,£)=0. 



Analogously to the treatment above, we introduce the form factors (i <E {d,N*}) as the eigenf unctions of the 
Lippmann-Schwinger kernel 



oo 

?n {£■) j=d,N* J 



(32) 



with an orthogonality condition analogous to (|29|). Then the separable form of the transition matrices is generated 
by the Hilbert- Schmidt expansion 

= J2w^(q,£)Q^(£)w^(q',£) , 0« {£) = f ^) 

In the spirit of the terminology, adopted for the separable approach, we may interpret the function Wi ;n as the form 
factor of the two-particle "bound state" j when the other two particles are in the "bound state" i. 

In Figs.|| through || we present the leading eigenvalues A n (£) and the real parts of r]n\£) and £n\£). Several 
comments are in order: 

1) The validity of separable expansion above is strongly limited to the energy region below the triton break-up 
threshold 3 H — > n + d, i.e., to energies £ < Ed, where Ed ~ —2.22 MeV denotes the deuteron binding energy. Above 
this region [£ > Ed), due to the singularities appearing in the propagators as well as in the potentials, the kernels of 
the equations become noncompact, and the Hilbert-Schmidt expansion loses its meaning. Therefore we shall restrict 
our consideration to the region below the triton break-up threshold. 



G 



2) The eigenvalue Ap goes through unity at the energy £ = £3 H sa — 13 MeV, which is essentially lower than the 
experimental value of the triton binding energy £^. p —8.5 MeV. This disagreement is a consequence of using the 
extremely simple Yamaguchi NN scattering matrix (||). Although this ansatz fits low-energy NN scattering, it has 
too long a tail into the high-momentum region and, therefore, predicts too much binding for the three-nucleon system. 

3) The modulus of the eigenvalues Aj and Aj+i are close to each other. Therefore, one can expect, that the leading 
terms in the expansion ( p2| ) will cancel each other to a large extent, which may lead to a relatively slow convergence 
rate of the sum. In order to demonstrate how well a finite sum (|22| ) may represent the exact amplitude, we show in 
Fig. H the ratio of the Schmidt norms 

Rn(£) = &, (34) 



of the operators 



TPQfi-V& n )y/i$>, (35) 



where V£[ n is given by the sum (|T^) containing only the first n terms. We see that the rate of convergence is not 
very impressive but appears to be sufficient for the practical calculation. Already in the region close to the two-body 
threshold (£ — e n ,d) the six leading terms in ([l9]) accumulate more than 95 % of the Schmidt norm of O. The same 
is valid for the rjNN case and especially for the {rjN) + (NN) channel as may be seen from Figs. [| and ||. 

The triton wave function was extracted from the pole of the 3iV scattering amplitude. In the present calculation 
we have restricted ourselves to the principal, spatially completely symmetric s-wave part, which in the usual notation 
(see e.g. J3l|]) reads 

* 3H = -r* s (Pi,P2,P3), (36) 

where £ a denotes the completely antisymmetric spin-isospin part, given by 

T^i=(x (0) C (1) -X (1) C (0) ). (37) 

The spin function x ( s = 0, 1) describes a three-nucleon spin=l/2 state which is obtained by coupling first the spins 
of two nucleons to a spin s and then by coupling s with the spin of the third nucleon to a total spin 1/2. The isospin 
function f '*) (t — 0, 1) is constructed in the same manner. 

The spatial part completely symmetrical with respect to the nucleon momenta, is given by 

y s {pi,P2,p 3 ) = ^{1 + Pi3+P2 3 )^(ki2,q 3 ). (38) 

Here, k\2 and q% are the usual Jacobi momenta for the tree (l+2)+3. The function i/j(k, q) can be expressed within 
our formalism in terms of the eigenfunctions ( p0| ) as follows 

i;(k,q) = ^(u^(k,q)-u^(k,q)), (39) 



with 



The normalization factor N is taken from the residue of the 37V scattering matrix (7) at £ — £3 H , and one finds 

-dXx 



N~ 



d£ 



F ■ (41) 



The elastic r]- 3 H scattering amplitude i^3jj(p) is then expressed in terms of the amplitudes X^l^, taken on the 



energy shell, 



F v »h(p) = (^i ( ;n (P- P> E ) + X i-il (P> P> E ) - 2X i]n (P- P> E )) ( 42 ) 



with p = A /2/i„3 H "E and ^,^3^ being the ?7- 3 H reduced mass. One should note that in the calculation of the amplitudes 
X\ entering (^), the driving terms Z\% and Zi\ (i = 2, 3) in (^|) and (^) were redefined in view of the approximation 
(|36|) (see Appendix B). 
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III. RESULTS AND DISCUSSION 



In Table || we present our results for the ry- 3 H scattering length obtained by keeping a finite number of terms in 
the Hilbert-Schmidt expansion of the amplitudes ([2^), ( p7f ) and (|3^). One can see, that the rate of convergence of 
a r) 3 H is approximately the same as that for the integral kernels discussed in the previous sections (see Figs. [$| f|and 
|). The choice n\ = 4, tt.2 = «3 = 6 provides rather satisfactory accuracy. Our result obtained with the ?yiV-input 
corresponding to a v N = (0.75 + ?0.27) fm [Q is 

a r) 3 H = (4.2 + i5.7)fm. (43) 

Firstly we note that the scattering length is rather large, which may be explained as a consequence of the virtual state 
(the real part of the scattering length is positive) generated by the strong attraction of the rj-3N interaction. The 
large magnitude of a ?) 3 H indicates that the corresponding pole of the scattering amplitude F v 3n(p) lies near the elastic 
scattering threshold. Turning to the region of positive energies, we present in Fig. || the total cross section for elastic 
ry- 3 H scattering. Here the virtual 7y- 3 H state manifests itself in a strong enhancement of cr(E) when E approaches the 
threshold value. It would appear natural that this effect is really observed, e.g., in ^-production in pd-collisions M. 
As already mentioned, we compare in Table [II our predictions on the ??- 3 H scattering length with those of 



using an r]N- interaction which reproduces the a^N scattering length of J^,|^] . One must note a rather strong difference 
of the results, especially for the real part of a^n- Apparently, the latter must be very sensitive to the position of the 
pole in the amplitude as it lies close to zero energy. Very rough agreement may be noticed with ||. But there is a 
tendency of our prediction to give a larger va lue o f 5Re a^3 H , which is probably due to the difference in the rjN range 



parameter r Q noted in the beginning of Sect. II B 



As to the question about the existence of a quasibound state, which an 77-meson can form with the lightest nuclei 
[p]-p|,p2"| , we see that within our model the r/-3N interaction is not strong enough in order to bind the ?7- 3 H system. 
This conclusion does not support the results of the FRA-model || or those of the optical model [§,[32| where an 77-3JV 
bound state appears already for rather modest values of the ijN scattering length. It should also be remembered that 
our calculations are based on the rank-one separable NN potential which is known to overestimate the attraction in 
the three-nucleon system. Therefore, we expect that the use of more refined NN models will most likely reduce the 
probability for binding the ^- 3 H-system. 



IV. CONCLUSION 



In the present paper, the four-body scattering formalism has been applied to study the r/-3N interaction in the 
energy region very close to the 77- 3 H elastic scattering threshold. The calculational scheme, which formally allows 
an exact solution, is based on the separable approximation of the appropriate integral kernels. The validity of this 
approach is justified by the fact that the driving NN and rjN interactions are governed mainly by the 5-matrix 
poles, lying in each case near the low-energy region. Therefore, the separable potentials, giving the correct structure 
of the amplitude close to the poles, are expected to provide an adequate approximation. In the present paper we 
have realized one of the possible schemes for solving the rj-3N four-body problem. To test its applicability, we have 
investigated the Hilbert-Schmidt procedure for constructing a separable representation of finite rank of the four-body 
kernels. The examination of their Schmidt norms shows that the rate of convergence by increasing the rank is not 
very fast, but satisfactory for practical purposes. Keeping 6 to 8 terms in the separable expansion, we have obtained 
a rather good precision. The same conclusion is valid for the calculation of the ?7- 3 H scattering length. At the same 
time we would like to point out that the present calculation suffers from the oversimplification of the iViV interaction, 
for which the rank-one Yamaguchi potential has been used. Of course, this shortcoming can be cured by using a 
more sophisticated separable approximation for the NN scattering matrix, requiring only more computational efforts. 
Thus, our results may be considered as a starting point for more realistic calculations. 



APPENDIX A 



In this appendix the transition amplitudes for r/-3N scattering are defined with due consideration of the Pauli 
principle for nuclcons. We begin with the two-body equations resulting from the quasiparticle approach to the four- 
body problem 0,0 

Z aTinn ,{f,p'\E)e 7n 4E)X 10]n , ln ,(p'\p\E)^ . (Al) 



■y^a. n" 
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Here the effective potentials 

Z a0 ., nn ,(p,p',E) =Y,<M£Mz)4 n ,(q',£') (A2) 

i 

are represented by the particle exchange diagrams (see e.g. Fig. ||). They include the propagator Ti(z) of the quasi- 
particle "i" depending on the two-body subenergy 

p 2 p' 2 p-p' 

2M a 2M P m k ' 

where M a denotes the reduced mass in the two-particle channel a. The form factors uf. n are generated by the 
separable expansion of the subamplitudes in the (1+3)- and (2+2)-partitions 

U%%q',£) =Y,u% n {m®an{£)ul n {q\8) . (A4) 

n 

In the case a = (1+3) the functions UQ(q,q',£) are the familiar Lovelace amplitudes appearing within the pure 
separable approach to the three-body scattering ^7j. For a = (2+2) these functions describe two independently 
propagating pairs of particles, each treated as a quasiparticle. In ( |A2| ) the energies £ and £' in the subsystems a and 
/3, respectively, are defined by 

2 /2 

£ = E- --— and £' = E - --— . (A5) 
2M„ 2Mb 



The momenta q and q' in (A2) are of course functions of p and p' 



q = p H p , q =p-\ p , (A6) 

rnfc m k 

whe re m k is the mass of the exchanged particle (or quasiparticle). After a partial wave decomposition, the equations 
( [Al| ) are reduced to a set of integral equations in one variable, being thus manageable for practical purposes. 

Turning to the r]-3N system we will firstly define the notation. The set of channels, wc arc interested in, is given 
by ([!]). In what follows, we do not consider the explicit structure of the equations, connected with the separable 



expansion of the basic subamplitudes (index "n" in (A4)) and drop also the spin-isospin indices. It is convenient to 
order the nucleons always cyclically, since in this case we have not to trace the sign, when changing from one state to 
another. Denoting the nucleons as particles 1, 2 and 3, we consider the following 7 states which contribute 

(i) one state in channel 1: |7y(123)), denoted as 1 , 

(ii) three states in channel 2: |l(ry23)), |2(?y31)), \3(r]12)), denoted as 2^ with i = 1,2,3 , 

(hi) three states in channel 3: |(1jj)(23)), | (297) (31)>, |(3fy)(12)), denoted as 3; with i = 1,2,3. 

Here we assume that the wave functions of the subsystems containing two and three nucleons are already antisym- 
metrized. For the present purpose we are interested in those amplitudes, which describe the transitions from channel 
1 to the channels 1, 2^ and 3;. They will be denoted respectively as X \, X ai with i — 1, 2, 3. Then using the separable 



approximation for the amplitudes (j22j), (27) and (|33|), the equations (Al) take the form 



3 3 

Xi = ^2 z - H -' v - + ^2 x - ■ (_) '' x '■ ' 

3 3 

.v, = z 2z . x + z, : :« .V; + 5^(1 - d ij )z 2i;2j e 2j x 2j + y ± -> H > A > . ( A7 ) 

3 = 1 3=1 
3 

= Z^-x + ^3 i; lOlXi + /.; j Hj A'_) , 

3 = 1 

where i — 1,2,3. The meaning of the driving terms Z CCi p j is explained schematically by the diagrammatical repre- 
sentation in Fig. ^. Their analytical expressions are given in Appendix B. The terms ^2,;3j an d Z-^^ij have different 
structures for i — j and for i ^ j and therefore are written down separately. 







The wave functions of different states belonging to the same channel may be obtained from each other by cyclic 
permutation of the nucleon coordinates. This fact allows one to obtain various relations between the transition 
amplitudes Z av ^.. For example, we have by cyclic permutation 

Z 2i;22 = (1(7723)1^2(7731)) = (2(7731)|^|3(77l2)) = Z 22 , 2s . (A8) 

Repeating the procedure, one hnds the general relation for i =/= j 

Z 2l ;2 j = Z 2l , 22 . (A9) 

In the same manner and by applying a combination of a cyclic permutation with a permutation within an antisym- 
metrized NN-st&te, one obtains for the transition 2^3 the general relations 

Z 2i;3i = Z 2l . M , (A10) 

Zl ::; Z 2 ::; , for i± J . (All) 

Then, taking into account the obvious identities 

e ai = e Q2 = e Q3 = e a , (a = 2, 3) (A12) 

and defining the properly antisymmetrized amplitudes in the channels 1 to 3 by 

i 3 

X? = Xi , Jf* = -^^X 041 for a -2,3, (A13) 
we arrive at the following set of equations 

xi = V3z 1;2l e 2 x 2 Q + V3z 1;3l e 3 x 3 a , 

A 2 a = V3Z 2i;1 + V3Z 2i;1 9 1 A 1 a + 2Z 2i;22 6 2 A 2 a + (2Z 2i;32 + Z 2i;3l )e 3 A 3 a , (A14) 
A 3 a = V3Z 3lll + V3Z 3l , 1 Q 1 X^ + (2Z 3l , 22 + Z 3u2l )0 2 X% . 

The last step we have to make is to change from the driving terms Z ai .p j to the terms, which couple the antisym- 
metrized states 

Zia = -^(77(123)1^1(7723)) = (77(123)|Z|1(7?23)) = Z 1;2l . (A15) 

Here P stands for the cyclic permutations of the nucleon labels according to 

\Pi{r)jk)) = \i(r)jk)) + \j{rjki)) + \k(r]ij)) . (A16) 



In a similar manner we obtain 



Z 22 = ~<Pl(»j23)|Z|Plfa23)) = 2Z 2l , 22 , 
Z 2 i = Z 2l -i , 

Z 13 = i=<77(123)|Z|P(l77)(23)) - Z 1;3l , (A17) 
Z 23 = i(Pl(7723)|Z|P(l77)(23)) = 2Z 2i;3a + Z 2l , 3l , 



^32 



2^3i;2 2 + -^3i:2i 



Substituting (A17) into (A14), we end up with the system 
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APPENDIX B 



In this appendix we list the explicit expressions for the driving terms, appearing in the separable-potential approach 
to the r/SN problem. As was mentioned already, we consider only s-wave orbitals in all subsystems. Therefore the spin 
algebra may be done independently. The corresponding spin-isospin coefficients S ss ' and A ss / are found by applying 
the usual angular momentum recoupling schemes (see e.g. p3| ). All driving terms have the same general functional 
form 

+1 

Z%L>(P*M = ?Y I ^Xq,e-^)T{E-f{p, P '))^{ q \£'-^)d{p-p'). (Bl) 



We list in Table IV for each driving term the corresponding assignments for the various symbols appearing in (Bl). 
where A ss < is defined by 

The coefficient 2 in the terms Z22 and Z| 3 stems from the identity of the nucleons, as described in Appendix A. The 
reduced masses are defined in (0), ([l5]) and (p4[). Furthermore, one should note that we have split Z23 into two terms 

Z23 — %23 + ^23 ■ (B3) 
The remaining driving terms are obtained from 

z&Z'M* E ) = zitS>Jp'>p> E ) > ( B5 ) 

Z&Z'M* E ) = 4tijP'>P> E ) ■ ( B6 ) 



In accordance with the approximation ( p6[ ) for the target wave function, we redefine the driving terms Z21 and 
Z31 in @ as well as the terms Z\% and Z\$ in (j^). Instead of them, we introduce effective potentials which are 
symmetrized over the singlet and triplet AW-states in the triton 

Z^ n ,(j>,p',E) Z { ° 2 %{p,p',E) = \ (zl 2 %,(p,pi,E) - z[%l,(p,p',E)) . (B7) 

— [ss ) — fss'J — (ss') 

The other terms Z 21 . nll , Z 13 . ln ,, and Z 31 . nll are defined by analogous relations. 
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TABLE I. 


Listing of parameters 


determining the separable 


parametrization of the AW and r/N scattering 


matrices. 


gf [MeV 1 / 2 ] 


[MeV" 1 / 2 ] 


I3 [fm" 1 ] j3i [fm- 1 ] 


g%, gj,, f3 v [fm- 1 ] fa [fm" 1 ] 


M [MeV] 


0.4076 


0.4863 


1.4488 1.4488 


2.10 1.04 6.5 4.5 


1661 



TABLE II. The scattering length a v 3 H as a function of the number n a of separable terms kept in the Hilbert-Schmidt 
expansion for the (1+3) and (2+2) amplitudes. The values of m, m and n-j, denote the number of terms kept in the separable 
expansion of the amplitudes (0) , (E?)) and (0) , respectively. 



m 


n-2 




a v 3 H [fm] 


2 


2 


2 


2.51 + i 1.60 


4 


2 


2 


2.54 + i 1.66 


4 


4 


2 


3.68 + i 3.46 


4 


4 


4 


3.99 + i 4.74 


6 


4 


4 


3.98 + i 4.75 


6 


6 


4 


4.16 + i 5.52 


6 


6 


6 


4.18 + i 5.67 


6 


8 


6 


4.19 + i5.69 
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TABLE III. Comparison of the predictions for the r\ 3 H-scattering lengths of the present model with the ones obtained in 
for the corresponding ryiV-scattering length. 



a vN [fm] 



S 3 H [& 



a v 3n [fm] (this work) 



0.57 + i 0.39 
0.29 + i 0.36 
0.27 + i 0.22 
0.55 + i 0.30 



1.32 + i 4.37 
0.58 + < 2.17 
0.41 +i 2.00 
-1.56 + i 3.00 



2.23 + i3.00 
0.97 + i 1.72 
0.69 + 2 0.67 
2.35 + i 1.68 



TABLE IV. Listing of symbols appearing in (Bl) 



n. 



r(ss') 



7 (ss') 



22;nn' 



^l(ss') 

z. 



23;nn / 
2(ss') 



23:nn' 



Sss< 

2A SS , 

Sss> 



» 



,,0) 



d;n 
..(«) 



P ' + I? 



-rl , 2 — 

P +sP 



01 _i Mjy^ — - 



p' + £tP 



p ' + X? 



'd 
T d 
TJV* 

rjv* 



2A d T 4Mj 



I PP 

~ 2M 



+ 



3p'" 



2n vN 4A/ N 

2 /2 

£ L L 

~ o r, _ _ ... i 



+ 



N 
P-p* ' 

P-P ' 



AT* 



2jS 



2fi riN 2p, d 



+ 



P*-p" 



P 



+ 



P I P-P 



d;n 

w 



, (*') 



P+2§tP' 



P+lTfp' 



P+ A/^+m^P 
P+^P' 



P+^P' 
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X, 



+ 



X, 



FIG. 1. Diagrammatic representation of the equations (pf) for the transition amplitudes X a of Tj-3N scattering. 
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FIG. 2. Diagrammatic representation of the potentials Z a p in the separable approximation (see Appendix B). 
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FIG. 3. Left panel: The leading eigenvalues of the Hilbert-Schmidt expansion for the three-nucleon scattering amplitude in 
the S — T — 1/2 state. Right panel: The ratio between the Schmidt norms of the kernels O and O n as defined by ( p^ ) and 
©■ 
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FIG. 4. The same as in Fig. § for ^AW-scattering in the (S; T) = (1; 0) state. 
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FIG. 5. The same as in Fig. | for the (<qN)+(NN) partition. 
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FIG. 6. Total cross section of elastic ?7- 3 H-scattering versus the cm. kinetic energy E. 
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